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SUMMARY 


An analytical model was developed for predicting the onset of supersonic stall bend- 
ing flutter in axial-flow compressors. The analysis is based on a modified two- 
dimensional, compressible, unsteady actuator disk theory. It is applied to a rotor 
blade row by considering a cascade of airfoils whose geometry and dynamic response 
coincide with those of a rotor blade element at 85 percent of the span height (measured 
from the hub). The rotor blades are assumed to be unshrouded (i.e., free standing) and 
to vibrate in their first flexural mode. The effects of shock waves and flow separation 
are included in the model through quasi-steady, empirical, rotor total-pressure-loss 
and deviation-angle correlations. 

The actuator disk model predicts the unsteady aerodynamic force acting on the cas- 
cade blading as a function of the steady flow field entering the cascade and the geometry 
and dynamic response of the cascade. Calculations show that the present model pre- 
dicts the existence of a bending flutter mode at supersonic inlet Mach numbers. This 
flutter mode is suppressed by Increasing the reduced frequency of the system or by re- 
ducing the steady-state aerodynamic loading on the cascade. The validity of the model 
for predicting flutter is demonstrated by correlating the measured flutter boundaries of 
several high-speed fan stages with their predicted boundaries. This correlation uses a 
level of damping for the blade row (i.e. , the log decrement of the rotor system) that is 
estimated from the experimental flutter data. The predicted flutter boundaries are in 
good agreement with the measured boundaries. 

These results show that the model can be used to estimate the relative stability be- 
tween operating points of a given rotor system. If, in addition, a measure of the mech- 
anical damping of the rotor system is available, the model can also be used to estimate 
the absolute stability at an operating point. 


INTRODUCTION 

Flutter can result in costly (both in time and money) overruns in turbofan-engine 
development programs. Solving the problem of flutter (at the engine development stage) 
may well mean major engine redesign and retesting. For this reason, engine manufac- 
turers and government agencies are currently supporting numerous research programs 
in an attempt to develop flutter prediction systems that can be used to design flutter-free 
engines. To date, these research programs have identified five regions (fig. 1) of the 


compressor performance map where flutter is generally encountered. 

Of the five regions shown in figure 1, the supersonic low-backpressure flutter re- 
gion has been the most thoroughly investigated analytically (refs. 1 to 3). In general, 
these analyses have considered the flow field through a cascade of two-dimensional air- 
foils undergoing simple harmonic pitching or plunging. The gas stream was assumed 
to be an inviscid, nonconducting, perfect gas. Shock waves that originated in the flow 
field were assumed to be weak so that supersonic small-disturbance theory could be 
used. Although these assumptions appear to oversimplify the flow conditions encoun- 
tered by a rotor at the onset of flutter, the flutter boundaries predicted by these analy- 
ses correlate well with experimental data. 

A recent analysis (ref. 4) has attempted to apply the supersonic, linearized, small- 
disturbance theory to higher backpressure operating conditions (i.e. , region IV of 
fig. 1) by including a finite-strength shock wave within the cascade passage. Results 
from this analysis show that the imsteady motion of the shock wave tends to induce bend- 
ing flutter. The existence of this flutter mode is documented in reference 5. 

The remaining operating region where supersonic flutter occurs (region V of fig. 1) 
lies close to the stall line of a stage . Analyses of this region have not appeared in the 
open literature. Literature on this subject (e.g. , refs. 5 to 7) generally presents ex- 
perimental data to document the extent of the flutter region. These data provide a very 
limited base from which an empirical correlation can be derived for predicting the onset 
of this flutter mode. 

The objective of the present analysis is to develop a model for predicting the onset 
of supersonic stall flutter. The calculated flutter boundaries of four high-speed fan 
stages that have experienced this type of flutter are compared with their measured 
boundaries to validate the analysis. 


BACKGROUND 

Supersonic stall flutter is generally encountered by rotors that do not have part- 
span dampers or tip shrouds. Experimental data reveal that the flutter mode is gener- 
ally the first flexural mode at a reduced frequency (based on tip relative velocity and 
tip semichord) of about 0.2. The vibratory pattern around the rotor tends to be very 
regular: All blades vibrate at the same frequency but are shifted in phase by a positive 
interblade phase angle of about 20° to 50° (ref. 7). This positive phase shift implies 
that the vibratory pattern is traveling around the wheel in the direction of rotation (i.e. , 
a forward-traveling wave) . 

A further characteristic of fans operating in flow regions associated with stall flut- 
ter is shown by the steady-state pressure distribution in figure 2. This pressure dis- 
tribution was produced from measurements taken across a rotor tip while the rotor was 
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operating near the flutter boundary. This figure clearly indicates a detached, leading- 
edge bow shock wave impinging on liie suction surface of the adjacent blade. The large 
extent of the compression region at the base of the shock wave seems to suggest a sepa- 
rated flow region that would increase in size as the fan operating point moved toward 
the supersonic stall bending flutter boundary. This separated flow would be accompa- 
nied by a decrease in shock strength. 

The present analysis develops a flutter model for the highly complicated flow field 
illustrated in figure 2 by using two-dimensional actuator disk theory. The effects of 
flow separation are included in the model through rotor-loss and deviation-angle corre- 
lations. For low-speed flows, actuator disk flutter models have been able to predict 
the onset of a single-degree-of-freedom bending flutter mode (e.g., refs. 8 and 9). The 
success of the models at low speeds suggests that a compressible actuator disk model 
might be capable of predicting the onset of bending flutter at supersonic speeds. This 
hypothesis has been confirmed by a number of calculations. These calculations, how- 
ever, did require as input the interblade phase angle at the onset of flutter. So that this 
requirement could be avoided, the actuator disk model was modified to allow for moder- 
ate values of interblade phase angle . This modification results in a flutter model that 
can be developed into a flutter prediction system. The validity of this flutter prediction 
system is demonstrated by comparing the predicted flutter boundaries of four high- 
speed fan stages with their measured boundaries. 


FORMULATION 

The present analysis for stall bending flutter of an isolated rotor models the rotor 
as a two-dimensional cascade of airfoils. The cascade is defined by the blade-element 
geometry on a cylindrical surface at a distance R from the axis of rotation. The flow 
field in the cascade plane is assumed to be two dimensional, compressible, and time 
dependent. Viscous forces are considered only within the blade channel. The unsteady 
flow variables associated with the rotor vibratory motion are assumed to be smaller 
than their steady-state counterparts. These variables at an instant in time are required 
to be periodic aroimd the wheel at a period equal to a fractional part of the circumferen- 
tial distance d = 27 tR. 

The motion of the airfoils in the cascade plane is restricted to simple harmonic 
plungi ng and edgewise motion at a cyclic frequency a? (fig. 3) . In addition the motion of 
each airfoil at an instant in time is assumed to be shifted from that of its neighbor by an 
interblade phase angle of a = 27m/N, where n is an integer and N is the number of 
blades in the rotor. In the present analysis the reduced-frequency parameter k|^ and 
the interblade phase angle are assumed to be small. Reduced frequency k|^ = wb/q_^. 
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where b is the semichord of an airfoil and q_^ is the magnitude of the relative inlet 
velocity in the cascade plane , ) 


Field Equations 

The governing equations for the flow field upstream and downstream of the cascade 
are written with respect to the relative coordinate system. The variables in these 
equations are nondimensionalized by the parameter q_^> the circumferential distance 
around the wheel d, the inlet steady-state fluid density p > and the inlet steady-state 
fluid temperature T_^. In addition the subscripts (oo, -<») signify upstream and down- 
stream flow variables. (The symbols used in this report are defined in appendix A.) 

The continuity equation for the upstream and downstream flow fields is 


9p- 


'±00 


at 


+ — (Uj 

dX, 


iboo^zboo 


) =0 


( 1 ) 


where p, U^, t, and are the nondimensional density, velocity components, and time 
spacial coordinates, respectively. (See fig. 3 for definition of coordinate system.) In 
addition the use of repeated indices denotes summation with respect to the repeated in- 
dex, The momentum equations for the flow fields are 


8U. . au. , 

Xy^oo ^ ttoo 


at 


8Xj 


1 

^±00 


ax. 


i = 1, 2 


( 2 ) 


where is the nondimensional pressure. The remaining field equations are the 

energy equation for an inviscid, nonconducting gas 


8s 


±00 


8s 


at 




±00 


ax. 


= 0 


( 3 ) 


and the equation of state 

s = L In P , In p, (4) 

y(y - l)Mf^ (y - l)Mf^ 

(The entropy was nondimensionalized with respect to T_^/q_^.) The variables appear- 
ing in equations (3) and (4) are the entropy sj the ratio of specific heats y; and the in- 
let steady-state relative Mach number M_ . 
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The preceding system of field equations is nonlinear. In this analysis they are cast 
into a system of linear equations with constant coefficients by expanding each equation 
to first order in terms of the amplitude of the perturbed variables . In the expansion 
process the zeroth-order terms are treated as constants . The resulting field equations 
are 




^±0 0 

at 


+ u. 




l,±oo 


'±00 _ 
+ A 
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axi 


±<x> 




axj 


(5) 


au. 


l,=boQ 


_ au, 


at 


+ u 




ax. 


1 

Pico 9Xj 


i = 1, 2 


( 6 ) 


9s; _ ar 

1,±00 


= 0 


at 
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(7) 


“ ±oo 


rsj 

A 


±00 


^ I 

±00 


( 8 ) 


m2 P±oo 

—OO ' — oo 

where the tilda denotes a time-dependent variable and the bar signifies a steady-state 
variable , 

The present analysis assumes that there are no sources of entropy or vorticity up- 
stream of the cascade. The upstream, unsteady velocity field must therefore be irro- 
tational and hence is equal to the gradient of a potential function , where 


a* 


Ui,-CO 


axj 


i = 1, 2 


(9) 


The unsteady velocity field dowmstream of the cascade is expressed as 

a$ 


u. 


1,00 


" + u? 


aXi 


1,00 


i = 1, 2 


( 10 ) 


where a3>^/aXj represents the irrotational component of the field and its rota- 

tational component. The source of the rotation field is the vorticity shed by the oscil- 
lating airfoils. The governing equation for this vorticity field is derived (ref. 10) by 
taking the curl of equation (6), which yields 


M+u. 
at 


_ai_ 
“ ax. 


( 11 ) 
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A solution of equation (11) is sou^t that is periodic in time and periodic along the cas- 
cade plane. This solution is given as 


f - 




( 12 ) 


where 






(13) 


U 


1, OO 


kg = 27m 


(14) 


(n is an integer) 


k = 


cud 




(15) 


and E is a complex constant. In reference 10 it is shown that a vorticity field has a 
vector potential associated with it. The relation between these two fields is given as 




(16) 




where V is the Laplacian operator. In addition it is shown that the components of the 
rotational velocity field are 


8X, 


(17) 


-§£- = {?« 
8X2 1 


(18) 


The particular solution of equation (17) with the right side given by equation (16) is 


* = E ^ikt 

'^oo 2 


R 


(19) 


Thus from ecjuations (17) and (18) the components of must be equal to 


Eiko 

Vy —e 

2^,2 

y 

^co ^ 


^^^3X2 ikt 


(20) 
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( 21 ) 


Solutions for the remaining flow variables are of the form 

^ ^ ^*^2^2 ikt 

P±co=%oo® ® ® 

S. =C, e”*le“‘2^2 3lfct 

±CX) zfcco 

S oo = 0 

r =Se'^le“'2^2^‘''* 

OO 

~ . qX ik2X2 ikt 

^±=0 = P±oo ® ® ® 

where C^_^, S, and p^^ are complex constants. Introducing these equations, 

along with equations (20) and (21), into the field equations yields 


( 22 ) 

( 23 ) 

( 24 ) 

( 25 ) 

( 26 ) 




p = -M^ (ik + Ut _ a + iU„ ^k„y e ^ (27) 

OO — OO y 1, — oo —oo Z / “cr« 
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a X, 


/ \ O' A- -iy 
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^ r. 

s = S e 
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-00 y 1,-00 -oo 2,-00 2 / — 




( 28 ) 

( 29 ) 

( 30 ) 

( 31 ) 

( 32 ) 
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( 33 ) 




P = -p (ik + U, a + U„ ik„'\c 

oo ^OO y 1,0000 2>00 2 iJ O 




where 


%«= 




]Aip 






1 - m|„(rq)u2 


1 - m|„(EQ)u2_^^ 


(34) 


and 


6 . = k + U„ . kn 

±oo 2 ,±oo 2 


RQ - 1 for < 0 




RQ for X- > 0 


(35) 


(36) 


In deriving these results it was assumed that the velocity field remained bounded at 
large distances from the cascade. 

The solutions for the unsteady flow variables are functions of four complex con- 
stants . These constants are evaluated in the next section. 


Boundary Conditions 

The constants C^_^, S, and E that appear in the solutions derived in the preceding 
section are evaluated here from four boundary conditions at the deforming cascade 
plane. These boundary conditions are derived by modeling the oscillating cascade as a 
deforming actuator disk. 

The first boundary condition requires the flow to be continuous across the deform- 
ing disk. The analytical form of this boundary condition is 

(l + p =(p at X =0 (37) 

y ooy n, -oo tx> ^coj n,oo 1 ' ' 

where represents the normal velocity components relative to the deforming 

disk. These velocity components to first order (i.e. , small cascade-plane deflection) 
are equal to 


8 
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(38) 
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OO 


8 X 
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fi. 

8 h- 

• Uo ^ — + 

84- 

8 t 

aXg 

8 X 2 


(39) 


where hj^ is the displacement of the disk from its mean position in the axial (X^) direc- 
tion (fig. 4) . The assumed analytical form for hj^ is 


ro 




e*“ 


(40) 


where h^^ is a real constant. Substituting equations (38) to (40) into equation (37) 
yields, to first order, 

U = p U at X- = 0 (41) 

-K30 ^OO OO X ' ' 

and 




— oo 





= "i(‘ 


ik + U„ ik„ - 

2,-00 2 


ik - 


P u. 

'^oo ( 


iko 


^ikt 


at 


Xj^ = 0 


(42) 


Introducing the expressions derived in the preceding section into equation (42) yields the 
algebraic equation 


“ll^-oo °^12^oo “l3^ “l4® 

where the constants *^12’ *^13’ ^14’ defined in appendix B. 

The next boundary condition requires the total enthalpy with respect to the deform- 
ing actuator disk to be locally conserved across the disk. The equivalent mathematical 
statement is 


y - 1 


T + i q^’ = — 1— T + - q£’ at X, = 0 

2 “ 2^ ^ 


(44) 


Rel 


where T^^ is the upstream and downstream gas temperature and q^_^ is the magni- 
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tude of the relative velocity approaching and leaving the disk. The expressions for 


are 


“—00 




u, + — 

8X, 




( 45 ) 



where is the local displacement of the cascade in the tangential (Xg) direction 
(fig. 4). The analytical form assumed for h.^ is 


rs^ 




‘‘^2^2 Ikt 
e e 


(47) 


where h 2 is the complex amplitude of the displacement. Combining equations (45) 
and (46) with equation (44) and expanding to first order yield 


+ +Up ) ^ + i(u? +uf ) atX^=0 

-1^2 2\ 1’-“° 2, -coy y.i~2 2\ 2,=o/ 1 


(48) 
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'84> 


84> 
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8X„ 8X, 



at X^ = 0 


(49) 


The local, total enthalpy change across the actuator disk as observed with respect 
to the Xj^,X 2 coordinate system is 


10 
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*^0, -oo 






1 1 
7 - 1 

-KXD 
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oO 
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+ U 


1 , 



_ 0$ 


+ u 


2, -oo 


ax„ 



at Xj^ = 0 


( 50 ) 


In deriving equation (50), variables higher than first order have been neglected. Com- 
bining equation (50) with equation (49) yields 



— oo 


Equation (51) shows that, if the steady aerod 5 mamic force is a finite quantity, the con- 
tribution of these forces to the total enthalpy rise across the disk is a first-order quan- 
tity. Combining equation (51) with equation (Cl5) of appendix C produces the second 
boundary condition 




8 t 




8 s 


dXr 
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0"^h, 


8 X 2 8 t 
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2 , 00 


O rsf 

8"^ho 


8X20t 


at Xj^ = 0 


(52) 


Substituting the solutions for the field variables into equation (52) yields 

“21^-00 °^2^oo ^ *^23® °'24® “ ^2 


(53) 


The coefficients in this equation are defined in appendix B. 

The third boundary condition requires that the local exit flow angle relative to the 
defoming actuator disk be specified. This angle is defined by the equation 


tan B 

^00 


n>j 



8 t 



0 t 


at Xj^ = 0 


(54) 
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Expanding the right and left sides of this equation to first order yields 


u„ 

tan B = 

U, 

1,00 


at = 0 


( 55 ) 


and 
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00 9 ^ 
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'a# 8h 
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where the local, perturbed exit flow angle /3 is defined as 


at Xj = 0 


(56) 


j 8 =8 -8 

^00 ^00 ^00 


(57) 


The present analysis assumes that the local, relative exit flow angle can be rela- 
ted to the local, relative inlet flow field (i.e., flow properties measured with respect 
to the deforming disk) through a quasi-steady deviation-angle correlation of the form 

6 = g^cv, a, 0'^ (58) 

where a is the local, relative incidence angle, is the local, relative inlet Mach 

number, ct is the local solidity of the deformed cascade, and 0' is the local staler 
angle measured from the normal to the leading-edge plane (fig. 4). The local, relative 
incidence angle is related to the local, relative inlet flow angle 8L00 I’^e^sured from 
the normal to the leading-edge plane by the equation 


a = 8 


1 

—00 



(59) 


The variable 6 in this equation is the angle between the chord line and the tangent 

mi 

to the camber line at the leading edge. The corresponding relation for the local devia- 
tion angle is 


2 

where is the local, relative exit flow angle measured from the normal to the 

trailing~edge plane and 6 is the angle between the chord line and the tangent to the 

m2 
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camber line at the trailing edge. If the leading-edge plane of the cascade is deforming 
in time, the local inlet and exit relative flow angles /3_^ and (3^ measured with re- 
spect to the Xj^-axis deviate from and j8^ by the angle tan”^ ^8hj^/9X2^, For 

small displacements, /3_^ and /3^ are 


B = S' + 

OO ^ — OO 


B = S' + 

^oo '^oo 




, ^ 

(61) 

8X2 


8h, 


^ 

(62) 

8X2 


Similarly the angle 0 between the chord line and the axial direction (fig. 4) is related 
to the local stagger angle S' by the equation 




0 = 0' + 


8X„ 


(63) 


Combining equations (58) to (63) produces the results 

a = B - 0 - 6 ^ 

m, 


8 = B^ - 0 - 

mr 


(64) 

(65) 


and 


^=0=^-^ 






( 66 ) 


In the present analysis the variables in equation (66) that are a function of the vibratory 
motion are B^o* ^_oo’ ^-oo’ Shj^/SXg. Expanding equation (66) to first order 

with respect to these variables yields 


^oo = ^ ^m„ + s( a. M_, a, 0' 




(67) 
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and 
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( 68 ) 


The local, relative inlet flow angle and the local, relative Mach number are de- 
fined by the equations 
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Expanding both these equations to first order yields 
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From figure 4, the local solidity of the deformed cascade is 
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Expanding equation (76) to first order yields 
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which for small interblade phase angles ^Ik^/N] « becomes 
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Introducing equations (72), (74), and (78) into the right side of equation (68) and 
replacing the left side with equation (56) produce the third boundary condition 
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Introducing the solutions for the field variables into this boundary condition transforms 
it into the equation 


“ 31^-00 ^ 32^00 “ 33 ^ “ 34 ® “ ^3 


(80) 


where the coefficients of the complex constants are defined in appendix B. 

The last boundary condition specifies the local entropy that is generated as the 
fluid passes through the cascade. The loss associated with the local entropy rise is de- 
fined in terms of a loss coefficient defined as 


X = 


pRel _ pRel 
0, -go 0,°o 



2, Rel 

'^—00 


(81) 


The variable iii this equation represents the local total-pressure loss 

measured relative to the deforming actuator disk. The entropy rise across the disk is 
related to this quantity by the equation of state 
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( 82 ) 


pRel _ pRel _ pRel 

0, "OO 0,OO 0, -OO 



Combining equation (81) with equation (82) and solving for yield 


s 

oo 




P-co<5- 


2, Rel 


0,-0 


(83) 


The present analysis assumes that the local loss coefficient x is related to the 
local inlet flow field, as observed with respect to the deforming disk, by a quasi-steady 
relation of the form 


X =X(«, M_^, ff, 0’) (84) 

The dependent variables in this equation are identical to those that appeared in equa- 
tion (58). Introducing equation (34) into equation (83) and expanding the resulting ex- 
pression to first order yield 


and 


s 

oo 9 





s = 





«Rel 

0.-°° 

^0,-oo 


1 

2 


2,RelN 



at X^=0 


(85) 


( 86 ) 


where is the local, perturbed inlet total pressure measured relative to the de- 

forming disk. This variable may be related to the local, perturbed inlet pressure 
and the local, perturbed inlet Mach number M_^ by expanding the definition of 


0, -oo 


= P-oo ^ 




M 


2,Reiy/<^-l> 


(87) 


to first order: 


17 






Substituting the solution for the flow variables into this last boundary condition re- 
duces it to 
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(92) 


C^A^C + 
41 “OO 


“ 42^00 


“43® 


^ 44 ® ~ ^4 


where the coefficients cn^j^, 0 :^ 2 * “ 43 * “ 44 ’ ®'4 defined in appendix B, 

The four boundary conditions applied in the present analysis have produced a sys- 
tem of nonlinear algebraic equations for the steady flow fields and four linear alge- 
braic equations for the complex constants E, and S. The steady flow equations 

are used to relate the steady downstream flow variables to the cascade geometry and 
tile steady inlet flow properties. These relations are incorporated into the algebraic 
equation for E, and S and render them independent of the downstream steady 

flow properties. The solutions for E, and S are combined with equations (19), 

(29), (30), and (31) and are used in the next section to calculate the unsteady aerody- 
namic force acting on the cascade. 


Aerodynamic Force 

The aerodynamic force acting on a deforming cascade is obtained by considering 
the flow field through a control volume that is fixed to a cascade passage (fig, 3). The 
present analysis assumes that the relative motion between neighboring airfoils, a mea- 
sure of which is the interblade phase angle, is small and that the reduced frequency of 
this motion is also small. The first assumption implies that the spatial variations of 
the flow variables across a cascade passage are small and thus can be neglected. The 
low-frequency assumption suggests that the rate of change of mass and momentum with- 
in the control volume is smaller than the net mass and momentum flux across the con- 
trol volume surface and thus can also be neglected in the present analysis. Based on 
these considerations the momentum equation for the control volume illustrated in fig- 
ure 3 can be written as 


F. = (p - p )n.T + - U,. ) (93) 

1 '^00 1 '1,00 1, —00' ' ' 

where is the force exerted by the airfoils on the control volume, t is the pitch of 
^ Rel 

the cascade, W is the mass flow exiting the control volume, and n^^ is the unit vec- 
tor normal to the exit plane of the control volume. Expanding equation (93) to first or- 
der yields 

= (p,. - p.„) V » - “1, -») 


and 
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( 95 ) 




F. = 
1 





HiT 


+ (p -p )E t+W^®Vu. - u. ^)+w(u. - u. ) 

Vt'oo ^-00/ 1 \ 1,00 1,-00/ \ 1,00 i,-oo/ 

The mass flow exiting the control volume is 


8h. 


W^® =p T U. +U. i n. 

^'oo \ 1,00 1,00 / 1 


(96) 


which, expanded to first order, becomes 


W = p tU. n. 

'^00 1,00 1 


(97) 


VRel 


8h. 


= p tU. n, + p tU. _n. + p r lU. i n. + p fUj _h. 

r^oo 1, 00 1 ^00 1, 00 1 ^00 1 1 af / 1 1 


8t 


(98) 


For the airfoil displacements under consideration in this analysis, the components of 
the unit vector n^ are 


rs/ 




(99) 


= -sin 

Expanding both of these equations to first order yields 




( 100 ) 


aoi) 
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Introducing these relations into equations (97) and (98) reduces them to 

W = p tU- 

'^OO 1,00 

and 


= p tU, + p tU, ^ + p t 

^oo 1,00 ^CO 1,00 ^oo 



_ 8h. 

p tU„ — ^ 


ao2) 


(103) 


(104) 


In addition if equations (101) and (102) are introduced into equation (95), the components 
of the unsteady aerodynamic force become 


rsj 




2 


(105) 

(106) 


The work done by the aerodynamic force on an airfoil over a cycle of motion is 
equal to 



(107) 


where /^e [ ] denotes the real part of the expression and the asterisk superscript signi- 
fies the complex conjugate of the variable. Inserting equations (40) and (47) into equa- 
tion (107) and evaluating the integral yield 


W. 

Aero 



(108) 
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where 9m [ ] represents the imaginary part of the expression wiihin the brackets in equa- 
tion (107). At the onset of flutter the aerodynamic work per cycle is equal to the mechani- 
cal energy dissipated over a cycle. If the mechanical dissipative force is proportional to 
the velocity of the airfoil, the mechanical energy dissipated over a cycle of motion is 



where the damping coefficient C is defined as 



2it d 


( 110 ) 


The variable 6 in this equation is the log decrement of the mechanical damping, and 
Myj is the mass per unit length of a rotor blade. This variable is defined as 

= Kp^TC ( 111 ) 

where K is a constant of proportionality, is the blade density, r is the average 
thickness of the blade, and c is its average chord length. Introducing this equation 
along with equation (110) into the integrated form of equation (109) yields 


W 


D 




hihi + h2h2j 


( 112 ) 


Equating equation (112) with equation (108) and solving for 6 establishes the minimum 
level of mechanical damping required for stability at an operating point as 


Onset 





Minj^ 9m(^FcJi 
2 _ 


2 ^ 2 ; 


(113) 


The symbol Min^^ denotes the minimum value of the function with respect to k 2 . If 
the available mechanical damping of a rotor blade exceeds this critical level, any small 
bending motion imparted to the blade will decay in time. Hence the system is stable 
and flutter will not occur. However, if is greater than the available mechani- 
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cal damping of a rotor blade, any small beading motion imparted to the blade will cause 
the blade to flutter. 

In the next section the present analysis is used to predict the flutter boundaries of 
four high-speed fan rotors , The results are compared with the measured flutter boun- 
daries to establish the validity of the analysis. 


RESULTS 

For tile flutter mode under study the amplitude of the motion increases monotoni- 
cally along the blade span from the node line. For blades that are rigidly fixed at their 
root (i.e. , mechanically constrained), the node line can be assumed to lie outboard of 
the blade platform (i.e. , outboard of the aerodynamic hub of the blade). The node line 
of the first flexural mode of the rotors to be investigated later in this section is approx- 
imately 20 percent of span height outboard of the platform. Thus, the vibratory motion 
of the blades in these rotors at any spanwLse section is 




h, = ? cos 


(114) 

ai5) 


where | is the amplitude of the motion at the section and e»(20%) ig the stagger an^e 
at 20 percent of span. 

The aerodynamic work per unit of span at a given radial location depends on the 
amplitude of the motion of the section and the relative dynamic pressure of the incom- 
ing streamli*ie to the section. Thus because both increase with distance from the hub, 
the total work done by the airstream on a blade is strongly influenced by the unsteady 
flow field surrounding the tip region, A simple calculation shows that this influence 
encompasses the outer 25 percent of the span. If the outer-casing boundary layer were 
assumed to influence the flow field over 5 percent of the tip region of a blade, a direct 
correlation might exist between the flutter boundary of a rotor and a two-dimensional 
cascade whose geometry and dynamic response coincide with that of the rotor at 85 per- 
cent of span (i.e., algebraic means of 75 and 95 percent of span). The results presented 
in the remainder of this section are based on this premise. 

The objective of the first series of calculations is to establish the influence of re- 
duced frequency k^ and interblade phase angle on the aeroelastic stability of a cascade 
of airfoils. For this study a normalized damping parameter defined as 
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6^ = - i ^ - sin 3 -ikt 

2 ^^0,Amb k^ ^ ' 


cos0*(2o%)j;J5;^3-''^2^2^-ikt 


is calculated as a function of k^ and the interblade phase angle. The parameter 
Pq in this equation is the inlet stagnation density measured in the absolute coor- 
dinate system. For greater than zero the airstream is supplying energy to the 

cascade. If the mechanical damping of the cascade system is zero, the operating 
condition would be unstable. For values of 6j^ less than zero the cascade is do- 
ing work on the airstream, and hence the system is stable. 

Solving equation (116) requires as input information the steady-state inlet flow 
properties, the geometry of the cascade, and the frequency of oscillation. A set of 
loss and deviation-angle correlations must also be specified. The present analysis 
makes use of the correlations derived in appendix C. The steady-state inlet flow 
conditions are derived from the blade-element data, measured at 85 percent of de- 
sign speed, of the second stage of the NASA 1450-ft/sec-tip-speed two-stage fan. 
(This rotor is analyzed in greater detail later in this section. This blade-element 
data set is reported in reference 6. The cascade geometry represents the geom- 
etry of the second-stage rotor element at 85 percent of span height from the hub. 

Figure 5 shows a plot of as a function of interblade phase angle and re- 

duced frequency kj^ for the NASA fan operating at the flutter boundary at 85 per- 
cent of design speed. This figure shows that the bending flutter mode is asso- 
ciated with a positive interblade phase angle (i.e., implies that the vibration pat- 
tern is traveling around the rotor in the direction of rotation) . The figure also 
shows that increasing the reduced frequency stabilizes the motion at a constant in- 
terblade phase angle. 

Figure 6 shows a plot of as a function of interblade phase angle and per- 

centage of design weight flow along the 85-percent-speed line. For a given inter- 
blade phase angle, decreasing the weight flow at constant wheel speed ^.e., in- 
creasing the steady-state aerodynamic loading) tends to destabilize the rotor. 

From the results shown in figures 5 and 6 the present analysis seems to pre- 
dict that increases monotonically with interblade phase angle over the range 

shown. Hence, unless some additional information is supplied that establishes the 
interblade phase angle at the onset of flutter, the analysis in its present form can- 
not be used to predict supersonic stall bending flutter. In attempting to solve this 
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problem it was observed that the variable kg (and hence the interblade phase 
angle) appeared in several of the terms in the analysis (e.g., eq. (77)) as a result 
of expanding the quantity 


/ ikg/N \ 
Nle ^ - l) 


(117) 


to first order in k„/N. It might therefore be expected that replacing k„ in the analy- 

^ /TVT \ ^ 


sis with the imaginary part of Nle 


/ ik„/N \ 
Nle ^ - ll 


9m 


N 




= N sin 


N 


(118) 


would give a more physically realistic behavior for large values of kg. This modifica- 
tion was introduced into all the expressions containing k„ except the exponential factor 
ik_Xg 

e ^ . The results obtained from this altered analysis are shown in figure 7 . For 

I ff| < 45°, the modified analysis yields results that are comparable to those obtained 
from the original formulation. In addition the modified analysis predicts a finite maxi- 
mum for 

The validity of the modified analysis as a correlation system for predicting the on- 
set of supersonic stall bending flutter can be established by showing that, for a given 
rotor, flutter will generally be observed whenever the predicted maximum value of 
exceeds the normalized structural damping of the rotor. For this demonstration, four 
high-speed fan stages were investigated. The level of structural ^.e., log decrement) 
damping (assumed to be a constant for each rotor) to be used in this study had to be es- 
timated from experimental flutter data. The first rotor to be analyzed in this study is 
the second stage of the NASA 1450-ft/sec two-stage fan. The performance map for this 
fan is shown in figure 8 . The second stage exhibits two zones of bending flutter, both of 
which are shown in figure 8 . The first flutter zone occurs at part speed and lies close 
to the fan stall line. The second zone occurs at operating speeds in excess of design 
speed. Because this second flutter zone appears to be associated with a choked passage 
flow, it will not be analyzed in this study. 

The maximum value of the normalized aerodynamic damping (denoted herein by 
the symbol nxax^ ^ function of weight flow and wheel speed for the second-stage 
fan is shown in figure 9. The operating point at 85 percent of design speed at which flut- 
ter occurred is denoted by the solid symbol. (Open symbols denote stable operating con- 
ditions.) If the level of mechanical damping in the system were equal to the aerody- 
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namic damping at this operating point, the altered analysis would predict that all the re- 
maining operating points are stable, including the stall points at 75 and 95 percent of de- 
sign speed. This result is supported by experimental results. An attempt was made to 
shift the flutter boundaries of this fan out of its stable operating range by closing the 
first-stage stator blades. This reduced the magnitude of the relative velocity entering 
the second stage and increased the relative incidence angle to the stage. Beducing the 
velocity raised the reduced frequency and, from previous results, is a stabilizing factor. 
However, increasing the incidence angle increased the aerodynamic loading on the rotor 
and from previous calculations, is a destabilizing factor. 

The computed results for max operating points (80 and 85 percent of 

design wheel speed) at which blade-element aerodynamic data were recorded are shown 
in figure 10. For a value of mechanical damping estimated from the results presented 
in figure 9, the analysis predicts that the modified design would experience bending flut- 
ter at these two operating points. Figure 11 shows the predicted flutter boundary for the 
modified fan superimposed on its performance map. The experimental flutter boundary 
is also indicated in this figure. Comparir^ the results presented in figure 11 with those 
presented in figure 9 shows that the analysis predicts a shift in the flutter boundary of 
the second stage as a result of the configuration change. The predicted shift of the 
boundary is consistent with the experimental data; however, the extent of the flutter re- 
gion is overestimated by the correlation. 

The present analysis can be used to study the effect of the steady, inviscid aerody- 
namic force on the supersonic stall flutter boundary. For example, reference 9 has 
shown that at low speeds the steady, inviscid aerodynamic force can induce bending 
flutter . In the present study, the aerodynamic damping was calculated as a function of 
interblade phase angle, and the effect of loss and deviation-angle variations was neglec- 
ted. This calculation was performed with, as input data, the rotor geometry and the 
aerociynamic inlet conditions at the predicted flutter point of the NASA two-stage fan. 

The results of this calculation are shown in figure 12, along with the results from fig- 
ure 7, which include the effects of the viscous forces. (The range over which the re- 
sults were calculated was chosen to be of sufficient width to include the interblade phase 
angles measured at the onset of flutter reported in ref. 8.) This figure clearly shows 
that the steady, inviscid forces in themselves are not the primary destabilizing force at 
high speeds . Rather it appears that supersonic stall bending flutter is induced by the 
viscous forces associated with a rapid rise in the tip blade-element loss coefficient. 

The present analysis was also correlated and verified with the flutter data obtained 
from fan C of the NASA Quiet Engine Program (ref. 6), the NASA 50-50 fan stage, and 
a prototype fen rotor developed in industry. The performance map for fan C (fig. 13) 
shows two flutter instability zones. The first zone is torsional flutter that occurs at 
part speed and lies close to the stall line. Just above this zone lies a zone of bending 
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flutter that appears to extend past 95 percent of design speed. Only the bending flutter 
mode of this rotor is analyzed in this report. The maximum value of the calculated, 
normalized damping as a function of wei^t flow and wheel speed for fan C is shown in 
figure 14. The operating points at which flutter occurred are indicated by the solid 
S 3 nnbols. (Open symbols denote measured stable operating conditions.) The level of 
mechanical damping for this rotor system is assumed to equal the minimum value of 
the aerodynamic damping computed at the operating point where flutter was encountered. 
This level of damping was used to calculate the stability boundary in figure 14. All op- 
erating points that lie above this boundary are assumed to be unstable, and those that 
lie below are assumed to be stable. The majority of the measured stable operating 
points lie below the estimated stability boundary. The results of transforming this 
boundary onto the performance map is shown in figure 15. The solid symbols in this 
figure denote measured flutter points. This figure shows that the present model yields 
results that are consistent with measured experimental data. 

The rotor of the NASA 50-50 fan stage encountered supersonic bending flutter at 
100 percent of design wheel speed at reduced inlet pressure (1/2 atm). At standard 
inlet conditions (1 atm) the flutter zone intersected the stall boundary at 90 percent of 
design wheel speed. Both boundaries are shown on the rotor performance map of fig- 
ure 16. The magnitude of the maximum normalized damping along the 100- and 
90-percent-speed lines for this rotor is shown in figure 17 . The experimental flutter 
point at reduced inlet conditions is indicated in the figure by the solid symbol. This 
point is assumed to define the mechanical damping level of the rotor system. Fig- 
ure 17 shows that the analysis predicts that all the stable operating points at reduced 
inlet conditions lie below the stability boundary. 

For a constant inlet-total-temperature operating condition the present analysis 
predicts that the aerodynamic damping acting on the NASA 50-50 fan stage will in- 
crease linearly with inlet total pressure. Thus if the mechanical damping of this rotor 
system were independent of the inlet total pressure, the predicted flutter zone on the 
performance map of this fan would increase as the inlet pressure was raised. In fig- 
ure 18 the predicted flutter boundary for the NASA 50-50 fan stage at an inlet pressure 
of 1 atmosphere is superimposed on its performance map. Also shown in this figure 
is the measured flutter boundary at that inlet pressure. Comparing these two sets of 
results with the measured flutter boundary at 1/2-atmosphere inlet pressure shown in 
figure 16 suggests that the present analysis overestimates the destabilizing influence 
of higher inlet pressure on this fan stage. 

A prototype fan rotor developed in industry was also run at two levels of inlet 
pressure. The flutter boundaries for this fan are shown on its performance map in 
figure 19, The maximum normalized damping for this fan is shown as a function of 
weight flow and wheel speed in figure 20 . The experimental flutter point for this fan 
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is denoted by the solid symbol . This point is used to define the stability boundary at 
atmospheric inlet conditions . Once again, the analysis predicts that all the measured 
stable operating points of this fan lie below this boundary. 

For a level of mechanical damping in the rotor system estimated from the results 
shown in figure 20, the analysis was used to establish the flutter boundary at a raised 
inlet pressure of 2 atmospheres. The predicted boundary and the measured flutter 
boundary at this higher pressure are shown on the performance map in figure 21. The 
measured flutter boundary intersects the stall line at 81 percent of design wheel speed; 
the predicted flutter boundary intersects the stall line at 71 percent of design wheel 
speed. As in the previous study, the analysis slightly overestimated the destabilizing 
influence of raised inlet pressure on the flutter boundary of a fan stage . 

The present analysis has been used with good success to predict the bending flut- 
ter boundaries of four transonic compressor rotors. These boundaries were calcu- 
lated from empirically derived estimates for the mechanical damping level in each 
rotor. The material damping level (i.e., log decrement) of a compressor rotor blade 
was reported in reference 4 to lie between 0.008 and 0.02. For the NASA two-stage 
fan the estimated damping level of 0.0126 lies within the material damping range re- 
ported in reference 4. This was expected since this rotor was designed for rig testing 
and thus was of fairly rigid and massive construction. For the remaining rotors the 
estimated levels of damping used in calculating their flutter boundaries were 0.0246 
for fan C, 0.0298 for the NASA 50-50 fan stage, and 0.0218 for the prototype fan 
stage. These levels of damping are greater than the level of material damping re- 
ported in reference 4 . The greater overall damping levels in these rotors are attribu- 
ted to their greater flexibility (i.e., these rotors were assembled from engine flight 
hardware) . This flexibility could allow the various components of the rotor assembly 
to rub against one another during a cycle of motion. This rubbing motion would re- 
sult in friction damping, which - added to the material damping of these rotors - would 
yield a total damping level greater than that in comparable rig rotor assemblies. Thus 
in using the present analysis as a correlation for predicting the onset of supersonic stall 
bending flutter of a compressor stage, a damping level of about 0.0125 should be used if 
the rotor assembly is fairly rigid. However, if the rotor assembly is fairly flexible, 
such that during a flutter cycle the components of the assembly might be in relative mo- 
tion, a damping level of about 0.025 should be assumed. 


CONCLUDING REMARKS 

In the present investigation an anal 3 d;ical model was developed to predict the onset 
of supersonic stall bending flutter in high-speed rotors. This model was based on the 
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assumption that there is a direct correlation between the onset of bending flutter in a 
rotor and the onset of bending flutter in a two-dimensional cascade whose geometry, 
dynamic response, and aerodynamic performance coincide with the rotor blade-element 
section at 85 percent of span height measured from the hub. The two-dimensional cas- 
cade of airfoils was considered as an actuator disk that was allowed to deform in time. 
The effects of shock waves and viscous forces were introduced into the model through 
rotor blade-element total -pressure-loss and deviation-angle correlations. The infor- 
mation input to the analysis is the steady aerodynamic inlet flow conditions to the rotor 
blade-element section at 85 percent of span height, the rotor geometry, and the rotor's 
dynamic response . In addition, aerodynamic design information must be specified, 
along with a value for the structural damping in the rotor assembly. 

To verify the analysis, computed flutter boundaries of four high-speed fan stages 
that experienced supersonic stall flutter were compared with their measured bounda- 
ries. These four rotors covered a range of design parameters typical of present and 
advanced fan designs. For this verification study the structural damping level had to 
be deduced from the experimental flutter data . 

The stability boundaries predicted by the analysis were in good agreement with the 
experimental results. Any difference that did exist showed the analysis to slightly over- 
estimate the extent of the flutter region. Although no specific reason for this apparent 
conservatism is obvious, several possible reasons are (1) the assumed structural 
damping model, (2) the use of rotor blade-element quasi-steady loss and deviation- 
angle correlations, (3) three-dimensional aeroelastic effects, and (4) the assumed form 
for the interblade-phase -an^e dependency used in the model. Finally the validation 
studies clearly showed that the primary destabilizing forces in supersonic stall bending 
flutter are the viscous forces associated with a rapid rise in the tip blade-element, 
total -pressure-loss coefficient. 

Additional experimental data are needed to improve the present correlation between 
the measured and predicted flutter boundaries of a fan stage and to broaden ttie general 
applicability of the model. In particular, experimental measurements are needed to 
establish the influence of rotor geometry, weight flow, and wheel speed on the flutter 
boundary and the interblade phase angle at the onset of flutter. In addition, a compre- 
hensive analytical and experimental program is needed to establish the mechanism and 
level of structural damping in a rotor assembly. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 9, 1978, 

510-55. 
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APPENDIX A 
SYMBOLS 

b airfoil semichord 

Cp specific heat at constant pressure 
complex constant 

c average blade chord and damping coefficient 
d circumferential distance around wheel 

E complex constant 

Fj aerod 5 Tiamic force 

h enthalpy 

hj^ axial displacement of an airfoil 
h 2 tangential displacement of an airfoil 

1 

K constant 

k reduced frequency based on circumferential distance, o)d/q_^ 

reduced frequency based on semichord, OJb/q_^ 

k 2 wave number of disturbance along cascade plane 

M Mach number 

m defined by eq. (C2) 

m^^ mass of a rotor blade per unit length 

N number of blades in rotor 

n integer 

n^ imit vector 

p pressure 

q magnitude of velocity field 

R radial distance 

RQ defined by eq. (36) 

R. complex constant 

S complex constant 


\ 
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^Aero 


entropy 
temperature 
total temperature 
time 

velocity components 
rotational velocity 

aerodynamic work per cycle of motion 

energy dissipated per cycle of motion by mechanical damping 
mass flow rate 
spacial coordinates 
incidence angle 

complex wave number in axial direction 
complex constants 
flow angle 

ratio of specific heats } axial wave number 
log decrement; deviation angle, eq. (37) 
metal angle at leading edge 

metal angle at trailing edge 

vorticity 

stagger an^e 

amplitude of blade motion 

density 

solidity 

cascade pitch or average thickness of a rotor blade 

velocity potential 

total-pressure-loss coefficient 

vector potential 

frequency of oscillation 
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Subscript: 


±oo upstream or downstream variable 

Superscripts : 

steady- state variable 
(~) time-dependent variable 

* complex conjugate 

Rel relative component 
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PERTURBATION COEFFICIENTS 
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APPENDIX C 


TOTAL-TEMPERATURE RISE ACROSS A DEFORMING ACTUATOR DISK 

In this appendix an equation is derived that relates the total-temperature rise 
across the deforming actuator disk to the upstream potential field and the downstream 
potential, vorticity, and entropy fields. The thermodynamic variables on each side of 
the disk satisfy the equation 


T ds = dh - ~ 


(Cl) 


where h is the enthalpy of the gas stream. Since equation (Cl) is a relation for the 
incremental change of state of the gas, it may also be written as 


1 T — 

p 8X. 8X. 8X. 


i = 1, 2 


(C2) 


Introducing equation (C2) into Euler's equation of motion 


1 , TT i - 1 8P 


8U. 8U. 

■+U. 

8t J 8X. p 8X. 


(C3) 


yields 


8U. 8U. 

8t ^ 8Xj 8X. 8X. 


i = 1, 2 


(C4) 


Substituting the relation 


8U. a /UiU.\ 

U. — L = -J-J. + U- , 

J 8X. 8X. \ 2 y J \ 8X 


8U. 8U. 

1 



(C5) 


into equation (B4) produces the result 
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For an ideal gas 
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is the total 


where Cp is the specific heat of the gas at constant pressure and Tq 
temperature. Introducing equation (C7) into equation (C6) yields 


8u. /8U. au/ 

— + vA—-- 

8t J 1 9Xj ax. 


T — - C ^ 

ax. P aXj 


(C8) 


This equation is often referred to as the Crooco-Vazonyi equation. The Xg compo- 
nent of equation (C8) for the flow field upstream of the actuator disk (in nondimensional 
form) reduces to 
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Similarly the relation for the flow field downstream of the deforming disk (in nondimen- 
sional form) is 
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Thus the difference in total temperature across the disk is equal to 
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1 8 


(C14) 



which, when expanded to first order, becomes 



(CIS) 
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APPENDIX D 


deviation-angle and loss correlation 

The deviation- angle correlation applied in the present analysis is based on a modi- 
fied Carter* s rule 

!> 

6 = 6^ + + r{a- aj^, M_^) (Dl) 

Va 

where cp is the camber angle of the blade cross section in radians and o!j^ is the in- 
cidence angle at the des:^ point in radians. The variable 6^ is the reference devia- 
tion angle for zero camber. This variable is assumed here to be a constant and is 
equal to the difference between the deviation angle at the desi^ point and the deviation 
ai^le calculated from Carter* s rule . The remaining variable m for a circular-arc 
airfoil is plotted as a function of stagger angle in reference 12 (p. 211). This graph is 
well approximated by the expression 

m = 0.216 + 0.05013 0* + 0.08617 0'^ (D2) 

where 0’ is ejqpressed in radians. 

In the neighborhood of the design point the deviation angle for a tip section of a fan 
rotor is generally observed to be independent of incidence angle and Mach number; 
thus 


r(“ - ^-cc) = 0 


Q!=o;t 


9r 

8Q! 


(a - ajj, M_^) = 0 


a=ol^ 
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9r 


8M 


(a - o-jj, M_^) = 0 


a=a-r 
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In addition the rate of change of the deviation angle with incidence away from the des^n 
point is generally observed to increase as the Mach number is increased. Based on 
these observations the expression 
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r = 


4,86 + - l) 

4.45 - 7,0(M_^ - l) 


\ 

r r 1 

/ \2 1 

u 

tan^Qfjj + o) + |l + 2 tan^ + dj 

(a - Q!jj) I 

1 

^ i 


(D6) 

was derived from a number of hi^- speed fan st^es tested at NASA Lewis, The data 
used in deriving this correlation were recorded at approximately 85 percent of span 
height (measured from the hub). The rotor-blade cross sections were multiple circu- 
lar arcs designed for high-speed operation (tip Mach number, ~1.4). The tip solidity 
of these rotors varied from 1. 3 to 1.7 . 

From equation (60) of the main text and the definition 

+ (D7) 

the quantity a - can be shown to be 

- “d " ^-co - ^=0 (1^8) 

Incorporating equation (D8) into equation (D6) and introducing the result, along with 
equation (D2), into equation (Dl) yield the deviation-angle correlation used in this 
analysis. 

The total-pressure-loss correlation is derived from a graph of experimental data 
presented in reference 11 (p. 248). This correlation is expressed as 


XyM?^ cos /3'^ 

f(D) 

4 a 



r/(r-i) 

- 1 


(DIO) 


where D represents the flow diffusion through a blade passage. For a two- 
dimensional compressible flow, D is defined as 


D = 1 + 


cos I3'_^ 
2cr 
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M 




cos 13' \ 
2(7 j 


(Dll) 
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An analytical expression for f(D) was derived by fitting a quadratic equation to the 
correlated experimental data of reference 11, measured at 85 percent of span hei^t. 
The result is 


f(D) = 0.03045 - 0.206 D + 0.48175 (D12) 

Introducing equations (Dll) and (D12) into equation (DIO) produces an expression for the 
loss coefficient x terms of inlet and exit flow angles, inlet and exit Mach numbers, 
and solidity. This expression, along with the deviation-angle correlation and the 
steady flow equations (eqs. (37), (44), (51), and (80)), forms a system of nonlinear 
algebraic equations whose solution yields the exit Mach number and exit flow angle as 
a function of inlet flow angle and Mach number, blade-section geometry, and design 
inlet and exit flow angles. From this solution the loss coefficient, as well as all the 
partial derivatives that appear in equations (73) and (86), can be calculated. 
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Weight flow 


Figure 1. - Compressor performance and stability map. 
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Figure 2, -Typical tip rotor static-pressure distribution near supersonic stall flutter boundary. 
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Figure 6. - Calculated, normalized damping as a function of interblade phase angle and weight flow for NASA two- 
stage fan at 85 percent of design wheel speed. 
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Figure 7. - Original and modified analyses for NASA two-stage fan at 85 percent of design wheel speed and 
75.68 percent of design weight flow. 
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Figure 8. - Performance map for NASA two-stage fan. 
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Figure 9. - Calculated, maximum, normalized aerodynamic damping as a function of weight flow and 
speed for NASA two-stage fan. 
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Figure 10. - Calculated, maximum, normalized aerodynamic damping as a function of configuration for 
modified NASA two-stage fan. 
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Calculated, normalized aerodynamic damping. 


3.4h- 


3.0 h- 


2.6 H 



i.sH 


1.4h- 


1.0 


30 


Wheel speed, 
percent of design 



40 50 60 70 80 90 100 110 

Weight flow, percent of design 


Figure 11. - Fan overall performance with first stator of NASA two-stage fan closed 5 degrees 
from design. 



Figure 12. - Effect of viscous forces on normalized aerodynamic damping for NASA two-stage fan at 85 percent of 
design wheel speed and 75. 68 percent of design weight flow. 
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Figure 15. - Fan C performance map showing stability boundary 
from figure 14. 
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Figure 18. - Predicted bending flutter boundary at atmospheric inlet conditions for NASA 
50-50 fan stage. 
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Figure 19. - Performance map for industrial fan stage. 


53 




Total-pressure ratio 



Weight flow, percent of design 

Figure 21. - Predicted bending flutter boundary at raised inlet conditions for prototype fan stage. 
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